Science Advances

Background: Pro-inflammatory and high-environmental-impact diets both threaten population and planetary health, but whether the two objectives align or conflict across countries is unresolved. We tested whether a supply-based dietary inflammatory index (sDII) is coupled to greenhouse-gas (GHG), land and freshwater footprints, and whether a nutrition-feasible reallocation can lower both simultaneously. Methods: From FAO Food Balance Sheets we built sDII (12 inflammatory-weighted components; construct validity r=0.9999) and five per-capita footprints using three independent life-cycle inventories for 182 national food-supply series. For each country, constrained optimisation reallocated 13 food-group supplies under isoenergetic, protein-preserving and food-group-bound constraints, minimising sDII and GHG jointly (Pareto frontier). Health burden was estimated via pooled relative-risk meta-analysis and 2023 World Bank population data. Results: sDII was only weakly associated with GHG (Spearman rho=0.14), land (rho=0.13) and freshwater (rho=0.28) in 2023. The balanced-Pareto reallocation lowered both sDII and GHG in 182/182 series (100% synergy): population-weighted delta sDII=-0.235, GHG -37.5%, land -49.3%, water -17.2%, i.e. 4.28 Gt CO2e/yr avoided. The associated reduction in metabolic-syndrome burden was directionally consistent but modest (~1.1% of the prevalent pool, ~2.84 million cases). Results were robust to three life-cycle inventories and three feasibility-bound regimes. Conclusions: Anti-inflammatory and low-carbon goals are decoupled rather than conflicting, and an isoenergetic, protein-preserving reallocation reconciles them in every country. Environmental gains are large and robust; health gains are directionally consistent but modest--triangulation, not a causal claim.

Species differ in longevity, physiology, and social organization, and these properties expose them to distinct endogenous and environmental mutagens. Mutational spectra generated by these processes can shape downstream molecular evolution, influencing synonymous nucleotide composition, codon usage, and even amino acid composition. Tracing this signal from life history to proteome through mutagenesis could reveal how mutational pressure interacts with the fitness landscape, including the direction of molecular change and the extent to which proteins remain functional while following mutational biases. Here, building on the recently identified age-associated mitochondrial A>G mutational signature in mammals, we test the universality of this signature and its downstream effects on genome and proteome evolution by comparing long-lived termites with short-lived non-termite cockroaches. We find that termite mtDNA exhibits a stronger A>G mutational signature than that of non-termite cockroaches, accompanied by coordinated shifts in synonymous nucleotide composition, codon usage, and amino acid composition. Our results show that ecological and life-history-associated mutational pressures can be transmitted through a hierarchy from mutational spectra to nucleotide composition and ultimately to proteome evolution. Mitochondrial genomes may therefore function not only as records of ancestry but also as molecular archives of the biological conditions under which species evolve.

Denisovans, originally identified from ancient genome from Denisova Cave in Altai, were a sister group to the Neanderthals and were once widely distributed across diverse terrains in the north and south of eastern Asia1-5. Genomic studies suggest that there were multiple events of interbreeding between modern humans (Homo sapiens) and Denisovans somewhere in Asia6. However, little is known about Denisovan living environments, diet, ecological niche, the timing of their disappearance, and the possible coexistence with modern humans in different regions. Here we report the radiocarbon age and stable isotopic signature of Penghu 3, a large Denisovan tibia from Penghu Channel, Taiwan7. The results showed that Penghu 3 dates to approximately 45,000 years ago, the time when modern humans were already widespread in southern parts of Asia. This Denisovan individual inhabited a C4-dominated ecosystem, open environments such as savannahs and floodplains, or a mixture of both, and consumed a high proportion of animal protein similar to some European Neanderthals8-10, with no clear evidence for the use of aquatic resources. These findings have implications for the behavioral flexibility, large body size7, and eventual disappearance of the Denisovans.

Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) drive stagewise membrane fusion by zippering into membrane-bridging four-helix bundles. Yet the conformations underlying successive fusion stages and the coupling of folding energy to bilayer remodeling remain unclear. Using optical tweezers, we measured the intermediates, energetics, kinetics, and force dependence of individual synaptic SNARE complexes assembled in cis on single membranes and in trans between apposed membranes. Membrane-anchored cis-SNAREs assembled through N-terminal and cooperative C-terminal/linker-domain transitions, whereas their transmembrane domains showed little intrinsic dimerization. Syntaxin retained membrane-dependent helical continuity through its linker domain before zippering was complete. PIP2 strengthened but slowed late zippering. In trans, membrane repulsion arrested single trans-SNARE complexes in a half-zippered state. G{beta}{gamma} further clamped this intermediate and inhibited late zippering; G-GDP, but not G-GTP{gamma}S, relieved the clamp, revealing a nucleotide-dependent mechanism for GPCR-mediated inhibition of neurotransmitter release. Modeling suggests that cooperative late zippering, syntaxin linker helicity, and concerted action of multiple SNAREs focus folding energy released over a long distance onto short-range membrane apposition. Thus, mechanically gated SNARE zippering is regulated by membrane forces, lipids, and regulatory proteins.

Osteosarcoma (OS) is the most common primary malignant bone tumor in children and adolescents. Despite aggressive multimodal therapy, including surgery and chemotherapy with methotrexate, doxorubicin, and cisplatin (MAP), outcomes for patients with relapsed or refractory disease remain poor, with no standardized second-line regimen. We developed a clinically calibrated, temporally resolved in vitro model to examine how resistance and collateral drug responses evolve during treatment with methotrexate, doxorubicin, and cisplatin (MAP). We modeled the evolution of chemotherapy resistance in vitro by exposing an OS cell line, MG63.3, to clinically relevant cycles of MAP therapy. MAP resistance and collateral drug responses were quantified over time across five independent evolutionary replicates, and three solvent-treated replicates served as controls. We complemented this phenotypic screening with transcriptomic profiling to extract predictive biomarkers of collateral drug response that could be used to personalize second-line treatment in patients with refractory OS. MAP exposure produced progressive resistance to doxorubicin and methotrexate, whereas cisplatin sensitivity remained comparatively stable. Repeated temporal screening against 12 additional drugs and combinations revealed three broad patterns of collateral response: progressive resistance, progressive sensitivity, and non-linear or stochastic change. Etoposide resistance emerged consistently, while sensitivity developed toward palifosfamide-etoposide and gemcitabine-docetaxel. Transcriptomic profiling showed dynamic, replicate-specific evolutionary trajectories rather than a single uniform resistance state. To demonstrate one approach for how our large, paired dataset could be used to translate these findings into potentially actionable clinical tools, we extracted predictive gene expression signatures of collateral drug response. By applying a dynamic, clinically inspired model of chemotherapy resistance in OS, we generated a detailed temporal mapping of collateral drug responses that highlights potential therapeutic windows and can be used inform selection of second-line agents. This approach and the resulting dataset can ultimately support personalized medicine strategies for patients with relapsed or refractory OS.

The evolutionary origin of the complete interferon (IFN) antiviral signalling cascade across chordates has long remained elusive, as fully functional IFN machinery had not been biochemically reconstituted in basal cephalochordates. Here we combine phylogenetics, AlphaFold3 structural prediction, multi-omic profiling and a full panel of in vitro and in vivo functional assays to characterize an intact primordial IFN network in amphioxus Branchiostoma japonicum, the extant sister group of all vertebrates. Amphioxus encodes a pair of undifferentiated BjTBK1/IKK{varepsilon} paralogs that form cytoplasmic heterodimers, whose conserved kinase-domain phosphorylation residues are mandatory for downstream IRF-driven IFN promoter activation. Despite minimal primary sequence homology with vertebrate IFNs, BjIFN1/2 possess compact -helical core folds that hint at potential structural analogy to vertebrate mucosal type III IFN-{lambda}. We identify two primitive class II cytokine receptors BjCRFB1 and BjCRFB2 that act as functional IFN sensors. Phylogenetic and tertiary structural validations confirm these undiversified ancestral receptor prototypes arose before vertebrate IFN receptor subfunctionalization, and their distinct temporal induction profiles upon viral challenge closely mirror the mucosal immune surveillance orchestrated by vertebrate IFN-{lambda}-IFNLR signalling. AlphaFold3 ligand-receptor docking identifies BjIFN2-BjCRFB2 as the optimal binding pair among all tested combinations, yet all modelled complexes yield low interface scores indicative of weak ancestral intermolecular interactions. Downstream signal transduction relies on two STAT paralogs, BjSTATa and BjSTATb, which jointly mediate transcriptional activation of core antiviral effector genes. The conserved GTPase BjMx serves as a central IFN-stimulated effector to suppress viral replication and maintain tissue immune homeostasis. Collectively, our structural and functional evidence demonstrates that the full hierarchical IFN signalling apparatus was fully assembled in basal cephalochordates prior to vertebrate radiation. The amphioxus IFN cascade preserves core ancestral molecular traits including structure-dependent ligand conservation and primitive low-affinity ligand-receptor interactions. This study resolves a long-standing evolutionary gap and provides definitive functional evidence demonstrating that IFN-mediated innate immunity constitutes an ancestral chordate trait, rather than a vertebrate-specific evolutionary innovation. It also puts forward the hypothesis that mucosal surveillance may represent the ancestral mode of chordate IFN defence.

Rib-cage morphology is a determinant of thoracic biomechanics, ventilation, and injury response, yet statistical shape models (SSMs) of the rib cage have relied on small cohorts (~100s of individuals) imaged by clinical computed tomography, which over-represents injury and disease. We constructed a surface-based SSM of the complete 24-rib cage from 26,275 standardised whole-body magnetic resonance imaging (MRI) scans of adults aged 19-74 years from the population-based German National Cohort (NAKO). Ribs were segmented with a deep-learning pipeline (a rib-extended SPINEPS model), reconstructed as per-rib surface meshes, and brought into dense vertex-wise correspondence by Gaussian-process morphable registration in Scalismo; the aligned ensemble was summarised by generalised Procrustes analysis and principal component analysis (PCA). Fourteen per-rib geometric descriptors provided a quantitative cross-walk between the abstract PCA modes and named shape features, and associations with sex, age, body size and composition (including body-fat percentage), and smoking exposure were estimated by multivariable regression with Benjamini-Hochberg false-discovery-rate control. Shape variation was strongly concentrated: 28 modes captured 95% of the total variance, and the first three alone accounted for 69.4% (PC1, 42.6%; PC2, 16.3%; PC3, 10.5%) and admitted consistent anatomical readings - a sexually dimorphic axis (PC1), a slender-versus-stout body-habitus contrast (PC2), and a free-rib-size axis at ribs 11-12 (PC3). The sexes were nearly fully separated along PC1 (Cohen's d = 2.52). Body mass and body-fat percentage were the dominant modifiable correlates of rib-cage shape, whereas the association with cumulative smoking exposure was comparatively small. The model is released as a population-representative geometric reference for benchmarking and morphing donor-derived finite-element human-body models and for further large-cohort shape analysis.

The membrane-bound DsrMK(JOP) complex is central to dissimilatory sulfur metabolism, including in sulfate-reducing microbes (SRM), a group that plays important roles in shaping planetary and human health. Despite this global importance, the mechanism of sulfide production and its links to energy conservation remain unclear. Here, we present high-resolution cryo-EM structures of DsrMKJOP from Archaeoglobus fulgidus, alone, with menadiol and with the sulfur-carrying substrate DsrC-trisulfide, complemented by physiological and biochemical studies. The results clarify how SRM control the reactivity of sulfur to selectively achieve sulfide production. While DsrC-trisulfide is highly stable in isolation, interaction with the DsrK subunit facilitates its hydrolytic activation, triggering a conformational change. This brings a perthiosulfenate sulfur intermediate into the catalytic pocket of DsrK for reduction at a single non-cubane [4Fe-4S] cluster, likely supported by a conserved non-ligating cysteine. DsrM harbors a structural quinone-binding site, but seems not to catalyze menaquinol oxidation, although this likely occurs in DsrMK complexes from different sulfur-metabolizing organisms. In DsrMKJOP, trisulfide reduction by DsrK is linked to quinol oxidation at DsrP, releasing protons to the periplasm to generate a proton-motive force.

The end-Permian and end-Triassic mass extinctions, driven by massive catastrophic volcanism and prolonged hypoxia, fundamentally reshaped life on Earth. However, the genomic impacts of these two biodiversity crises on living organisms remain largely unknown. Here, we performed a genome-wide screening to identify accelerated evolved regions in the ancestral lineage of mammals that survived both extinction events. Nearly all (20/21) of these accelerated regions were located in protein-coding sequences, and 81% (17/21) were found to be associated with lung function. We further extended our analysis to three additional vertebrate lineages that experienced either one or both of the mass extinctions. Similar genomic signatures, involving accelerated evolution of lung-related genes, were also observed in these non-mammalian lineages. Collectively, these findings suggest that adaptation of lung-related genes to prolonged hypoxia may have occurred during the two mass extinction events, potentially driving a second evolutionary stage of the lung following the vertebrate transition to land.

HP6/Umbrea, a rapidly evolving Drosophila (Heterochromatin Protein 1) HP1-family paralog, has well-documented sequence and regulatory evolution but an under-studied molecular function. In this manuscript, we hypothesize that HP6/Umbrea acts as an HP1-recruited plasticizer, providing support for this model using coarse-grained molecular-dynamics simulations of HP1a condensates. Retaining only the dimerizing chromoshadow domain (CSD), HP6/Umbrea notably lacks independent chromatin-binding capacity but binds HP1a directly, co-localizing with it in vivo. We report that when covalently tethered to an HP1a carrier, HP6/Umbrea partitions into HP1a condensates ~6-fold more strongly than when free, supporting HP1a-mediated recruitment as its entry route. Once incorporated, HP6/Umbrea leaves the phase-separation threshold, interfacial tension, and host partitioning statistically unchanged, but monotonically lowers dense-phase density. These observations are consistent with a spacer function rather than generic loss of cohesion. Importantly, unchanged short-time internal mobility suggests a packing effect, predicting increased permeability to large transcriptional machinery, potentially resulting in a position effect-variegation (PEV)-like modulation of heterochromatic silencing. Finally, comparative sequence analysis shows the C-terminal tail is a recently originated, purifying-selection-constrained innovation, which is consistent with an evolved function in this region. In sum, our simulations suggest a mechanistic basis for how HP6/Umbrea may have evolved as a condensate plasticizer and thus potentially act as a rheostat for leaky transcription.

Soil microalgae endure harsh terrestrial stressors, such as intense light. Eustigmatophytes are an independent evolutionary branch within stramenopiles and occupy diverse aquatic and terrestrial environments, but the structural organization of their photosynthetic apparatus remains poorly understood. Here, we determined the cryo-electron microscopy structure of a photosystem I-light-harvesting complex I (PSI-LHCI) supercomplex bound with ferredoxin-NADP+ oxidoreductase (FNR) from the terrestrial eustigmatophyte Vischeria stellata at 2.44 [A] resolution. The supercomplex contains a monomeric PSI core associated with only three LHCI subunits, representing the smallest PSI-LHCI reported among structurally characterized red-lineage PSI complexes composed of violaxanthin-Chl a proteins (VCPs). The three VCPIs with distinct structure features and arrangements form a compact belt along the PsaL-PsaI-PsaM side of PSI. The structure also resolves a 43-residue N-terminal segment of FNR (FNR-N) bound to the PSI stromal surface, which is stabilized by both a eustigmatophyte-conserved insertion in PsaL and the N-terminal region of PsaD. In contrast, the catalytic region of FNR was not resolved, suggesting conformational flexibility. Computational simulations indicate potential excitation-energy-transfer pathways connecting the three VCPI subunits to the PSI core and highlight lineage-specific pigments that maintain energetic connectivity within the exceptionally compact antenna. Our analysis further reveals conservation of FNR tethering despite pronounced diversification of antenna size and organization. These findings uncover a modular evolutionary principle in which PSI acceptor-side organization is retained while the light-harvesting antenna is extensively remodeled, providing a framework for understanding the diversification of photosynthetic energy conversion across ecological transitions.

The insect functional homolog of the mammalian vagus nerve has remained elusive. Here, we identify the SIFamide (SIFa) peptidergic circuit as this homolog in Drosophila melanogaster. SIFa neurons project the longest axons from the brain to the hindgut ampulla, forming synapses with peripheral SIFa-receptor (SIFaR) neurons originating in the abdominal ganglion. This circuit is functionally plastic: sexual experience enhances synaptic strength at the hindgut, increasing excretion while reducing mating duration to to favor energy conservation and metabolic recovery. Mechanistically, this shift involves octopaminergic/tyraminergic gating and SIFa-mediated modulation of Adipokinetic Hormone (AKH) from the corpora cardiaca. We further uncover a feedback loop in which SIFaR/Allatostatin-A (AstA) neurons project back to the brain, signaling hindgut physiological status to SIFa neurons via AstA-R1. This hindgut projection is conserved across diverse insect orders. Thus, the SIFa-SIFaR circuit constitutes an evolutionarily conserved "insect vagus nerve" that integrates metabolic state and reproductive history to optimize survival strategies.

Classical synaptic transmission assumes that amino acid neurotransmitters and monoamines are stored in distinct vesicle populations. Here, we provide structural evidence for a dual-transmitter organization in which histamine and GABA are spatially partitioned into a core and rim within individual synaptic vesicles. Using an engineered glutaraldehyde-NaBH4; immuno-electron microscopy platform, quantitative analysis of 2,195 vesicles revealed a previously unrecognized nanoscale architecture, with histamine-associated signal concentrated within the vesicular core (95.1%), while GABA is organized toward the peripheral rim. Light microscopy further demonstrated extensive histamine-GABA correspondence across central and peripheral tissues, including sympathetic ganglia and adrenal medulla. This intravesicular segregation challenges the conventional separation of amino acid and monoamine storage into distinct vesicle classes and reveals that chemically distinct transmitters can occupy organized domains within a single vesicular lumen. This architecture may provide a structural basis for the distinct physiological modes of GABAergic and histaminergic signaling, linking nanoscale vesicular organization to their established differences in temporal action.

Visual neuroscience requires precise spatiotemporal projection of optical stimulation onto the retina, especially in experimental mouse models. However, in vivo patterned stimulation in mice is profoundly hindered by the extreme optical power and severe anatomical aberrations of the eye. Consequently, visual stimulation relies mainly on unverifiable, open-loop approximations that often lack spatial precision. Here, we introduce a closed-loop, spatially modulated stimulation platform that overcomes these barriers. By integrating a digital micromirror device (DMD) with electronically tunable lenses (ETLs) and a real-time, fundus camera-guided focus optimization module, we directly verify the location of patterned stimuli on the retina while dynamically correcting for chromatic and geometric defocus. This platform delivers quantitatively verified static and dynamic patterned stimuli to the living retina with lateral resolutions as fine as 6.7 {micro}m. Guided by ray-tracing optical analysis, our work establishes a technological foundation that enables highly reproducible, cellular-scale interrogations of the visual pathway.

Stable nitrogen isotopes of amino acids are widely used to reconstruct trophic position. Most applications rely on only two amino acids despite routinely measuring many others. Here, comparative amino acid {delta}15N values from 88 mammal species reveal that threonine records a physiological dimension beyond trophic position. Adding threonine to the canonical glutamate-phenylalanine framework reveals ecological differentiation obscured by broad dietary categories and opposite isotopic relationships between herbivores and secondary consumers. A mechanistic model links this variation to preferential intestinal utilization of threonine for mucin synthesis and predicts experimentally testable patterns of isotope partitioning. These findings show that amino acid {delta}15N values encode complementary dimensions of organismal ecology, expanding amino acid isotope analysis beyond trophic reconstruction to reveal physiological and ecological variation associated with dietary specialization.

Photosystem II (PSII) is in eukaryotic phototrophs is generally considered to operate within a more restricted spectral range than photosystem I (PSI), in which long-wavelength chlorophylls are a well-established feature of the peripheral antenna. Whether eukaryotic PSII can acquire comparable far-red-associated properties through lineage-specific antenna diversification has remained unclear. Here we present a 3.09 [A] cryo-electron microscopy structure of the C2S2M2L2 PSII supercomplex from Euglena gracilis, a euglenophyte species harbouring a secondary plastid and unusual light-harvesting system. We show that the euglenophyte-specific antenna protein LhcE9 occupies the position corresponding to canonical Lhcb5, but in a markedly different orientation that creates a distinct interface with the PSII core, particularly with CP43. Combined structural, spectroscopic, mutagenesis and proteomic analyses support LhcE9 as the stably bound PSII antenna subunit most closely associated with the far-red state in the supercomplex. Excitation-energy-transfer calculations further indicate two fast lineage-specific antenna-to-core routes mediated by LhcE9 and PsbX. Together, these findings reveal an unexpected mode of PSII antenna diversification and provide a structural framework for far-red-associated light harvesting in PSII.

Accurate accounting of aquatic methane emissions is critical for climate change mitigation, yet current global budgets overlook a key driver: the elevation-regulated atmospheric pressure. Here, we present the first large-scale investigation of methane ebullition across 164 shallow waters spanning elevations from sea level to 4886 meters. We demonstrate that ebullition rate increases with elevation-over four times higher at >3000 m a.s.l. than at sea level-due to two synergistic, pressure-dependent physical mechanisms: a degas effect (enhanced bubble formation) and a trigger effect (facilitated bubble ascent). Independent theoretical prediction of the combined effects shows near-perfect agreement with the empirical elevation trend, quantitatively confirming that these physical mechanisms are the primary drivers of enhanced ebullition at high elevations. Our findings reveal that mountain aquatic ecosystems represent unaccounted methane hotspots that have been systematically underestimated in global inventories.We therefore call for urgent integration of these ecosystems into IPCC assessments and targeted mountain mitigation and sustainable management strategies.

Sensory neurons have been increasingly recognized as vital contributors to deep tissue function. However, how these specialized neurons contribute to salivary gland function remains largely undefined. Here, we uncover a role for trigeminal somatosensory afferents in salivary gland perception and function using in situ-based classification, in vivo calcium imaging, behavioral assays, and targeted ablation. Retrograde labeling from the submandibular gland complex revealed substantial direct innervation from trigeminal neurons. Further categorization confirmed that Trpv1+ sensory neurons provided dense innervation of the Whartons ducts. TRPV1 agonist ductal infusion directly activated gland complex-associated neurons in the trigeminal ganglia and evoked a robust pain phenotype. Targeted Trpv1+ ablation disrupted Whartons ducts structure and dramatically reduced stimulated saliva volume. Our work provides the first evidence that Trpv1+ sensory neurons maintain salivary architecture and are necessary for stimulated saliva production, revealing a vital interoceptive role for direct trigeminal innervation in submandibular gland health.

Denisovans are an extinct archaic Homo group whose lineage diverged from the Neanderthal lineage approximately 550,000 years ago and were widely distributed across eastern Asia until [~]45,000 years ago1-7. Their morphological features are known directly from the existing cranio- dental and phalangeal remains1,2,8-12. However, the body size and postcranial morphology of the Denisovans remain largely unknown. We here report that hominin femoral and tibial fossils recovered from the Penghu Channel, Taiwan, are Denisovans in their proteomic profiles. Morphologically, these specimens are among the largest leg bones known in Pleistocene Homo. They exhibit generally archaic features, but also show some modern human-like morphology, including a strong femoral pilaster. Our findings demonstrate that the Denisovan population at the northern circle had larger body size than earlier Homo erectus as well as Late Pleistocene Homo sapiens in eastern Asia. This challenges the generally held expectation that Pleistocene Homo followed Bergmanns rule that anticipates latitudinal decline of body size, and suggests that the large Denisovan brain resulted from their large body size at least partly. The strong pilaster developed in the Penghu femur suggests some behavioral similarities between the Denisovans and the Upper Palaeolithic modern humans and/or gene flow from the latter to the former.

Previous molecular archaeological research has made the surprising discovery that almost all studied European walrus finds dating between the 11th and 14th centuries CE seemed to be traded from the Norse colony of Greenland. But why has almost no walrus ivory from the Barents Sea region been detected in medieval Europe, despite its closer proximity and being mentioned in historical accounts? Here, we show that the full spatial extent of medieval Atlantic walrus hunting has been underestimated in previous analyses that targeted modified walrus skulls (rostra) with which pairs of attached tusks were traded, because destructive sampling of ivory artefacts is not always practicable. By analysing workshop offcuts from medieval Sigtuna, Sweden, we directly compare the geographical sources of walrus ivory with rostra, and show that isolated tusks were traded in different ways. We first resolve the genome-wide trans-Atlantic walrus population structure, discovering significant genome-wide nuclear differentiation between western and eastern Atlantic walrus populations. By employing spatially diagnostic nuclear SNPs (n = 144,986), we then use low-coverage sequence data to classify historical walrus rostrum and tusk specimens to either western or eastern Atlantic origin, and thus overcome previous limitations when provenancing walruses based on mitochondrial DNA alone. We find that all medieval walrus rostra are assigned exclusively towards the western Atlantic. In contrast, half of the medieval ivory specimens from Sigtuna are assigned to eastern Atlantic sources including Iceland and the Barents Sea region. Moreover, the objects of eastern origin precede those from Norse Greenland in time, suggesting sequential exploitation. Our observations resolve a discrepancy between earlier molecular inference and historical evidence and imply a significantly broader extent and impact of medieval ecological globalisation in the Arctic.